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Abstract 

The hadronic cascade description developed in an earlier paper is extended to 
the response of an idealized fine-sampling hadron calorimeter. Calorimeter response 
is largely determined by the transfer of energy E w o from the hadronic to the elec- 
tromagnetic sector via 7r° production. Fluctuations in this quantity produce the 
"constant term" in hadron calorimeter resolution. The increase of its fractional 
mean, f2o = (E^o) /E, with increasing incident energy E causes the energy de- 
pendence of the 7r/e ratio in a noncompensating calorimeter. The mean hadronic 
energy fraction, ff L = l — f® , was shown to scale very nearly as a power law in 
E: /° = (E/Eq)™- 1 , where E « 1 GeV for pions, and m ps 0.83. It follows that 
7r/e = l — (1 — h/e)(E/ Eq) 171-1 , where electromagnetic and hadronic energy de- 
posits are detected with efficiencies e and h, respectively. If the mean fraction of 
/° which is deposited as nuclear gamma rays is f®, then the expression becomes 
7r/e = 1 — (1 — h'/e)(l — f®)(E/Eo) m ~ . Fluctuations in these quantities, along 
with sampling fluctuations, are incorporated to give an overall understanding of 
resolution, which is different from the usual treatments in interesting ways. The 
conceptual framework is also extended to the response to jets and the difference 
between it and p response. 
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1 Introduction 



In Paper I[l] we developed a conceptual basis for understanding the division 
between hadronic and electromagnetic (actually tt°) energy deposition in a 
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Fig. 1. Energy flow in a hadronic cascade. A fraction f n o (with energy-dependent 
mean /° ) is transferred to the electromagnetic sector through ir° production in 
repeated hadronic inelastic collisions. The ir° and hadronic energy deposits after 
the division are separately stochastic, and so must be treated as parallel statistical 
processes. Each produces a visible signal, whose sum E V1S is sampled. 

contained hadronic cascade. 1 The model "calorimeter" was a very large iron 
or lead cylinder, with no energy leakage except via muons, neutrinos, and 
front-surface albedo losses. Extensive Monte Carlo simulations gave results in 
good agreement with test-beam measurements. The relevant conclusions of 
the paper were that: 

(1) All significant hadronic energy deposition is by low-energy particles 

1 GeV), whose energy and species distribution in a given medium is 
independent of the energy or species of the incident hadron. (Hadronic en- 
ergy was denned as all energy not carried away by ir° decay photons.) The 
existence of this "universal low-energy hadron (and nuclear gamma ray) 
spectrum" makes it possible to define an energy-independent efficiency h 
for the conversion of this energy into a visible signal in a fine-sampling 
calorimeter. 

(2) In each high-energy collision of the hadronic cascade, a significant fraction 
(typically 1/4) of the energy is lost from further hadronic activity via 
7T° production. A sequence of high-energy hadronic collisions bleeds off a 
larger and larger fraction of the energy as the incident energy E increases. 
The net fraction transferred to the ir° sector in a given cascade is f n o, and 
the mean tt° fraction is f® - 2 This one-way flow is illustrated in Fig. 1. 3 

(3) In particular, the mean fraction of the energy /° = 1 — /° in the hadronic 
sector scales very nearly as a power of the incident energy, 

f h (E) = (E/Eor- 1 , (1) 



Most of the content of Paper I was first presented at the 1989 Workshop on 
Calorimetry for the Superconducting Super Collider [2]. 

2 Throughout this paper, the superscript indicates the mean of a stochastic vari- 
able, e.g. f® = (/no). In most of the literature the quantity without the superscript 
indicates the mean. 

3 Wigmans points out that the actual number of 7r°'s produced is quite small[3]. 
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where m ~ 0.83 (with some mild absorber Z dependence) and Eq 
1 GeV for pions and ~ 2.6 GeV for protons (again, with some Z de- 
pendence). Physically, m is related to the mean number of secondaries 
and the mean energy fraction going into 7r°'s in any given collision in 
the cascade, and Eq is the energy at which multiple pion production be- 
comes significant. Both must be determined by experiment for a given 
calorimeter. 

(4) It was predicted that a calorimeter would have a different response to a 
proton than to a pion. 

The observations pertain equally well to a homogeneous or fine-sampling 
calorimeter, and have significant implications for its response and resolution. 
"Fine-sampling" means that absorber and sensor elements are thin compared 
to both the em radiation length and the neutron interaction length. It has 
the same structure throughout: no separate front em compartment or rear 
catcher. It can be an inorganic crystal calorimeter, a uranium/liquid argon 
calorimeter, or a lead/scintillator- fiber calorimeter. 

The power-law approximation given in Eq. (1) is just that, for reasons discussed 
in Paper I. It seems to work well over the energy range of available test-beam 
data, about 10 GeV to 375 GeV, and it has the required asymptotic properties: 
It is everywhere positive, and /° — > {n/e — > 1) as E — > oo. The physical 
assumptions it is based upon become less dependable at very high energies and 
are not valid at energies below the threshold for multiple pion production. 

As far as possible, results in this paper are obtained without recourse to the 
power- law approximation for /°, in order to obtain more general results than 
those relying on this more approximate form. 

Most of the results reported in this paper can be found in the Proceedings 
of various conferences and workshops[2,4-9]. The Monte Carlo results used 
in these papers are often based on now-superseded versions of hadronic cas- 
cade simulation codes[10], the oldest being FLUKA86. In particular, nuclear 
gamma rays were not included, so that the em deposit is exclusively via ir° 
production. Since these versions many improvements in the codes have been 
made, e.g., improvements in FLUKA by Ferrari and Salafll], especially in 
the nuclear physics modeling. The failings of the old code are apparent in 
Fig. 2(b), for example, where the points fall below the 45° line because of 
unscored hadronic energy. A large fraction of the unscored energy is evidently 
that of nuclear gamma rays. On the other hand, ir° energy deposition was 
very well described [12] and can be trusted. In Paper I we reported simula- 
tions with MARS10, HETC, and FLUKA, which, though based on different 
high-energy interaction models, were in excellent agreement. Since in this pa- 
per I depend only upon the high-energy division between the 7r° and hadronic 
sectors, calculations based on the older code have not been repeated. 
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In Sec. 3, I distinguish between em energy deposit by ir° decay photons and 
by nuclear de-excitation gamma rays. A fraction f n o of the energy is deposited 
via 7T° decay, and a fraction f^f^ by nuclear gamma rays within the acceptance 
gate, where fh = 1 — f^o. The total em deposit is E n o + = E(f w a + fhf-y)- 

Recent developments are incorporated, some of which were predicted or dis- 
cussed in Paper I. These include Cherenkov readout [13], which is sensitive 
only to em (ir° and nuclear gamma) energy deposition, and observation of the 
n/p response difference [14]. 

Central to the paper is the discussion of resolution, where conditional prob- 
ability distribution functions (p.d.f.'s) are combined to account for parallel, 
independent stochastic processes. 

Hadron calorimetry is a well-traveled road, explored in hundreds, if not thou- 
sands, of papers over several decades. The object here is to present a broad- 
brush treatment of hadronic cascades in a simplified generic calorimeter, in 
hopes that a somewhat nonstandard approach can contribute to our physical 
understanding of a real calorimeter. Real calorimeters, with front em compart- 
ments, rear catchers, leakage, crack corrections, jet finding algorithms, and a 
myriad of other features, are described in dozens of test-beam study results, as 
well as in published studies of compensation, the role of neutrons, and other 
matters. These are discussed in detail in Wigmans' book[3] and review[15], 
the review by Leroy and Rancoita[16], and in their many citations. None of 
these practical problems are discussed here. 



2 Albedo and f^o 

The 7T° fraction f n o = E n o/E increases with energy, but at any given energy 
it is subject to large fluctuations. FLUKA simulations of the 7r°/hadronic en- 
ergy division are shown in Fig. 2. The model absorber consisted of a large lead 
cylinder (50 cm radius, 250 cm long) in which the first 25 cm (about 1.5 inter- 
action lengths) was treated as a separate region. In Fig 2(a) no distinction is 
made between the regions, while in 2(b) interaction of the incident pion was 
not permitted in the front section, but energy deposited there is included. It 
acted as a catcher for back-scattered interaction debris. The distribution about 
the ideal E^ = 1 — E n o shows less scatter in 2(b) because front-face, or albedo, 
losses are included. Most albedo loss comes from backward or backscattered 
products of the first collision; when the first interaction occurs deep in the de- 
tector there is essentially no albedo loss. Runs at 50 GeV with and without an 
"albedo catcher" show an average difference in deposited energy is 0.43 GeV, 
or 0.8%. Out of 1000 cascades 50% lost less than 0.2 GeV, and 3.4% lost 
more than 2 GeV. In the simulations the amount of lost albedo energy rises 



4 



30 



25 



> 

O 20 



15 



10 



25 



> 

O 20 



S 15 



c 

1 io 



5 - 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 

30 GeV %~ on lead 
(a) Albedo energy lost 



I I I I I I I I I I I I I I I I I I I I I I 



(b) Albedo energy included 



5 10 15 20 25 
Electromagnetic energy (GeV) 



30 



Fig. 2. Calculated energy deposit distributions for 30 GeV negative pions incident 
on a lead "calorimeter." In case (a) backscattered energy is lost; in (b) it is retained. 

only slowly with increasing incident energy, as might be expected. While these 
losses are not totally negligible, I omit them from resolution considerations in 
Sec, 7 because (a) the distribution is sharply peaked at near-zero loss, and (b) 
the losses are small, particularly at higher energies. 

For reasons discussed in the introduction, the points shown in Fig. 2 scatter 
below the 45° line because older versions of FLUKA did not account for all of 
the hadronic energy deposit, even in the absence of albedo losses. Presumably 
most of this downward scatter (about 15% in the worst case) is the result of 
the program's failure to tally nuclear gamma rays, most of which come from 
de-excitation following slow neutron capture by nuclei. According to Ferrari 
and Sala[17], these might account for nearly 10% (Fe) or 20% (Pb) of the 
7T° + 7 fraction, or 5%-10% of the total energy deposit. These contributions 
scale with the hadronic fraction, not the ir° fraction. Given that the hadronic 
fraction is underestimated by this fraction in this simulation, it is better to 
take the hadronic fraction as 



/h = 1 - />, or E h = E-E n 



(2) 



5 




Fig. 3. FLUKA simulations for negative pions incident on a lead "calorimeter." (a) 
Distribution of the ir° energy deposit for 20 and 200 GeV incident pions, and (b) 
energy dependence of the mean ((f n o) = f® ), standard deviation (a n o) of the f n o 
distribution, and standard deviation relative to the mean (cr^o/fo). The dotted line 
is discussed in the text. 

As the number of Monte Carlo events in the sample increases, the (E^^E^o) 
distribution projected onto the E n o axis approaches the marginal distribu- 
tion Il(/ 7r o), the p.d.f of the 7r° energy fraction. Two (unnormalized) FLUKA- 
generated examples of f n o = E-^o/E distributions are shown in Fig. 3(a). The 
mean and standard deviations are shown in Fig. 3(b). 

The fractional mean /° moves slowly to the right with increasing energy, and 
can be represented by /° = 1 — (E/Eq)" 1 " 1 . As it does so, the rms width of 
the distribution decreases only slowly (presumably because of increasing 7r° 
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multiplicity), and is well-represented by 4 

a n o = 12.5% x (£/100 GeV)~ - 076 . (3) 

There is no physical basis for this functional form, except that it remains 
positive as E — > oo. As the incident energy becomes very large, the distribution 
"crowds" the right limit, and the variance should approach zero. The fits to 
real data (Fig. 14 and Tbl. 3) yield somewhat different values for the multiplier 
and exponent, which in any case should vary from case to case. 

Alternatively, one might have chosen to express the width as a fraction of the 
mean rather than as a fraction of the total incident energy, or as a n o / /° rather 
than a n o. A power-law fit to the Monte Carlo data in this form is indicated by 
the dashed line in Fig. 3(b). The strong energy dependence just reflects the 
energy dependence of /° , and the mildness of the energy dependence of cr n o 
is obscured. 

The dotted line in Fig. 3(b) is a fit by Acosta et al. to SPACAL data for 
Crma(/i)//i in the range 10-150 GeV, where "1" refers to the central tower [18]. 
The authors show that the f\ distribution is a good representation of the f n o 
distribution. The fit is given as 0.435 — 0.052 ln(E). SPACAL was a lead/scin- 
tillator fiber calorimeter, while our model is a solid lead cylinder. The test 
beam events contained contributions from nuclear gamma rays. Even so, the 
difference between the dotted line and the FLUKA-simulation dashed line is 
difficult to understand. 

The dimensionless "coefficient of skewness," 71 = /i 3 /crj (where /13 is the third 
moment about the mean), is constant to within the Monte Carlo statistics with 
a value near 0.6. There are no significant higher moments within the sensitivity 
of the simulations. One might expect to see the skewness change sign; the tail 
should move from the right to the left side of the most probable value at very 
high energies. This transition has not yet been observed. 

It is instructive to examine a continuous distribution with similar properties: 
the Beta distribution f(x;p, q) = x 9 ~ 1 (l — x) s_1 / B(p, q), where the normaliz- 
ing constant B(p, q) is the Beta function[19]. It may be thought of as "a con- 
tinuous version of the binomial distribution:" It is defined only for < x < 1, 
is zero at both limits, and the peak position, variance, and other properties 
depend on p and q. For p > 2 and q > 2 it has zero derivative at both limits. 
If the mean is less than 0.5 the distribution is skewed to the right, as is the 
case for U(f n o); for larger means the distribution is skewed to the left. The 
skewness of II(/ 7r o) remains positive for means much greater than 0.5, but, as 
explained above, might be expected to change sign as the mean approaches 

4 A slightly better fit is obtained with 0.126 - 0.0099 ln(E/100 GeV). 
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Fig. 4. The projection of "events" onto the diagonal line at 6 = tan -1 (h/e) gives 
the energy distribution for a calorimeter in the absence of any other fluctuations. 
Events must fall along the thick solid segments. One can imagine the projection as 
increases from 0° (h/e = 0, large dotted histogram through the undercompensated 
case shown in (a)) to 45° (h/e = 1, full compensation) where it approaches a 
(5-function. For 9 > 45° (h/e > 1, shown in (b)), the calorimeter is overcompensated, 
and the skewness changes sign. 

unity. 

Different or better cascade simulations would be expected to produce distri- 
butions with somewhat different shapes and different moments. What is of 
consequence here is that a function H(E n o) exists which describes the energy 
distribution of the 7T°'s for a given primary energy E; no significant conclusions 
in this paper depend upon the details. 

The contribution of U(f n o) to the calorimeter resolution can be understood by 
a geometrical construction. Figure 4 shows the same MC "events" as Fig. 2(b), 
but with the lost hadronic energy restored as per Eq. (2) (except for some 
vertical scatter retained for clarity). The observed energy distribution in the 
absence of sampling fluctuations is the projection of this distribution onto 
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a diagonal line at 9 = tan~ l (h/e). The limits of the projected distribution 
[f n o = and f n o = 1) are shown by the arrows. All of the events thus project 
onto the solid segment of the line, with length | cos 9 — sin 9\ = (1 — h/e) cos 9. 
The sampled n(/ 7r o) distribution in Fig. 3(b) is replotted as the histogram 
along the f n o axis. A point at f n o = 1 projects to the end of the solid segment, 
so the energy scale along this axis is foreshortened by cos 9. The length of the 
solid line segment, rescale by 1/ cos^, is (1 — h/e). The fractional standard 
deviation of U(f n o), a n o, also scales as (1 — h/e). It thus contributes (l—h/e)a n o 
(in quadrature) to the calorimeter resolution. 

For h/e = 1, 9 = 45° and the distribution becomes a 5-function. For h/e > 1 
(overcompensation) the distribution "flips," with the tail on the low-energy 
side, since the 7r°-rich events in the high-energy tail of Il(/ 7r o) now contribute 
less energy to the cascade than do the hadrons. 

Experimental verification of this situation is at least strongly suggested by the 
WA 78 results obtained with a uranium-scintillator plate calorimeter [20]. The 
bulk of the energy was deposited in the upstream "Section A," which in some 
configurations was Fe-scintillator and in others U-scintillator. Unweighted en- 
ergy distributions are shown in the paper's Fig. 3 for the Fe-scintillator case 
(Fe25) and in Fig. 8 for the most uniform U-scintillator case (U15). While other 
resolution effects broaden the distributions, the distributions are nonetheless 
skewed to the right for Fe25 and skewed to the left for U15. The dotplots 
in their Figs. 4 and 9 show uncorrected .E^t vs A max , the maximum energy 
deposited in one of the scintillator sheets. Large deposits indicate large em 
shower activity, and hence 7r°-rich events. The mean of the distribution slopes 
upward in the Fe25 case and downward in the U15 case, again providing evi- 
dence for the inversion of the distribution. 



3 ir/e 

An electromagnetic shower initiated by an electron or 7r°-decay photons pro- 
duces a visible signal (potentially observable via ionization or Cherenkov light) 
in a calorimeter with efficiency e. Most of the ionization is by electrons and 
positrons with energies below the critical energy, of order 10 MeV (21.8 MeV 
for iron, 7.0 MeV for uranium). The response, here temporarily called "e," is 
usually linear in the incident energy E, and so serves to calibrate the energy 
scale: 

"e"=e£ (4) 

As shown in Paper I, the visible signal produced by hadron interactions also 
comes predominately from low-energy ionizing particles whose spectra and rel- 
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ative abundance are independent of the incident hadron energy. Many mecha- 
nisms are at play, including endothermic nuclear spallation. Neutrons play an 
especially significant role [21]. These mechanisms are exhaustively treated in 
the literature; for example, in Refs. [3,15-17]. The sum of all the hadronic en- 
ergy deposit mechanisms (excluding showers by 7r° decay photons) produces 
an observable signal with efficiency h. In most cases h/e < 1. For a mean 
hadronic fraction = 1 — /° , 

V= e&E + hfiE 

= eE[l-(l-h/e)f h ] . (5) 

In the case of an an incident pion, the response relative to an electron is 
*/e = 1 - (1 - h/e) ft . (6) 

Specializing to our power-law form for /°, 

vr/e = 1 - (1 - h/e){E/E Q ) m - x = 1 - aE m ~ x , (7) 



where, as above, m ~ 0.82 to 0.86. (ir/e is only defined for an ensemble of 
events, so it is implicitly a mean value.) Since the physics leading to the power 
law involves a multistep cascade, it is not expected to be dependable below 
5-10 GeV. Only m and a = (1 - h/e)E^ m can be obtained from fits to data, 
at least in a single-readout calorimeter. For incident pions (not protons) a 
range of Eq's near 0.7-1.0 GeV fit almost as well because Eq is raised to a 
small power. The ratio h/e cannot be obtained from a measurement ofn/e as 
a function of energy without other information or some assumption about Eq. 

I emphasize again that the power-law representation is not empirical, but 
follows from an induction argument. It has the correct asymptotic limit, since 
f® — > 1 as E — > oo. For a 10 19 eV proton-induced air shower, for example, 
/°o ~ 0.98, in accord with the usual cosmic ray expectation and observation 
that nearly all the energy deposit at very high energies is electromagnetic. 5 
The expected behavior of 7i/e is shown in Fig. 5. 

Representative fits of test-beam results to Eq. (7) are shown in Fig. 6. Solid 
curved are least-squares fits with both m and a allowed to vary, while dashed 
curves are fits with m is constrained to its nominal value from Paper I. Given 
the wide range of experimental data which have been fitted to test the power 



5 This is an illustrative example only, because there is no expectation that m will 
remain even relatively constant over such a large energy range. In addition, as much 
as 10% of the energy is carried by muons and neutrinos from meson decay [22]. 
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Fig. 5. Energy dependence of ir/e expected for several values of h/e with the as- 
sumption that Eq = 1 GeV. For almost all calorimeters, h/e < 1. The value for a 
given combination (U/scint, etc.) depends on the actual configuration. The lower 
dotted line, for h/e <C 1, should be applicable to a calorimeter with Cherenkov 
readout. The power-law description is not expected to be dependable below about 
10 GeV, but nonetheless seems adequate at 5 GeV. 



Table 1 



Power law fits to a variety of ir/e measurements. The ZEUS uranium/scintillator[23] 
and DO U/LAr[24] calorimeters are so close to compensating as to be uninteresting 
in this context. 



Calorimeter 


m 


a 


x 2 


h/e* 


Expected h/e 


SPACAL[25] 


0.788 


0.164 


9.2 


0.836 


o.sss-o.sgst 


(lead / scint-fiber ) 


0.830$ 


0.141 


14.0 


0.859 


0.853-0.895+ 


CDF end-plug had cal[26] 


0.865 


0.244 


2.7 


0.756 


0.667 t 


(50 mm Fe/3 mm scint) 


0.816$ 


0.286 


14.1 


0.714 


0.667 t 


Copper /quartz-fiber [13] 


0.833 


0.753 


2.6 


0.247§ 




(QFCAL) 


0.816$ 


0.814 


3.8 


0.238 




U/scint (WA 78) [20] 


0.85* 


-0.555" 




1.555* 





* Assuming E = 1 GeV. 

f Paper I, Table 1. (The calorimeters have only approximately the same structure.) 

$ Dashed curves in Fig. 6: m held fixed at the value given by the fitted line in 

Paper I, Fig. 12. Error in m from this work is ±0.01 to ±0.015. 

§ Akchurin et al. report e/h ~ 5 [13]. Virtually all of the hadronic Cherenkov 

signal can be accounted for as coming from relativistic pions[27]. 

(t Ref. [20] gives two data points, shown in Fig. 6. Here I assume m = 0.85 and 

adjust a for a best fit. 
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Fig. 6. Fits to test-beam results for a lead/scintillator-fiber)[25], for the CDF up- 
grade end-cap hadron calorimeter (50 mm iron/3 mm scintillator sheets) [26] , for 
a copper /quartz- fiber test calorimeter [13], and for the WA 78 uranium/scintillator 
calorimeter [20]. Fit parameters are given in Table 1. 

law, one suspects that occasional disparate results (e.g., the low value of m 
for the CDF end-plug calorimeter) indicate data reduction problems. 

Although the energy fraction carried by the nuclear gamma rays scales as 
the hadronic fraction, it is detected with efficiency e (nearly). 6 Let / 7 be 
the fraction of the hadronic energy deposited by nuclear gamma rays within 
the electronic gate time. Its mean, /°, is independent of both incident hadron 
energy and species via the "universal spectrum" concept developed in Paper I. 
Thus fh f-y of the incident energy is detected with efficiency e, and the total em 
fraction / em is f n o + fhfy The remaining fh (1 — / T ) is detected with redefined 
hadronic detection efficiency h! . Then 

vr/e = (/°o + f h f") + (h'/e)f° h (l - /°) 



6 The Compton electrons are low-energy to begin with, and a significant fraction of 
their energy deposit occurs by ionization after they drop below Cherenkov threshold. 
They are thus detected with lower efficiency than are high-energy electromagnetic 
cascades. This nicety is ignored in the present discussion. 
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l-(l-h'/e)(l-$)fi 
(l_fc'/e)(l-/j) {E/E T~ l 

l-aE m - 1 , (8) 

where, as elsewhere, the superscript zero indicates the mean. 

The important point here is that the power-law description given by Eq. 7 
is recovered, even though part of the electromagnetic signal tracks with the 
hadronic sector. It is immaterial whether we use (1 — h/e) or (1 — h'/e)(l — 
/°), except that we should remember that h contains a nuclear gamma-ray 
component. 

Energy deposit by nuclear gamma rays can account for 10%-20% of the total 
energy deposit in materials such as iron, and it can be even higher in high- 
Z materials [17]. Since most of the nuclear de-excitations are the result of 
slow neutron capture, nearly all of the gamma rays are emitted on a time 
scale of hundreds of ns (see Fig. 3.22 in Ref. [3]). Electronic gate widths for 
calorimeter signals are made as short as possible, given the arrival time spread 
of information in such large devices, and this turns out to be close to 100 ns. 
This means that most of the nuclear gamma energy deposit happens outside 
of the sampling time. There are a few gammas from faster neutron captures, 
but by and large most of the gamma signal is lost. Thus / 7 as defined above 
is in the range of a few percent. While the above development is of interest 
in showing that the power-law scaling does not need modification, it is of not 
much practical importance. 



4 n/p 

We observed in Paper I that /° is larger for an incident charged pion than 
for an incident proton (or neutron). This is a consequence of the fact that a 
leading hadron, carrying a large fraction of the energy, is likely to have the 
same quark number as the incident hadron. If the collision is instigated by a 
charged pion, there is high probability that the leading hadron is a 7r°, but for 
an incident proton or neutron the leading hadron is most likely a baryon. 

Via the "universal spectrum" concept we expect that the hadronic response 
of the calorimeter is the same for protons as for pions, except for a scale 
factor: h/e is the same for both cases, and the mean hadronic fraction ratio 
/°- / fp < 1 is independent of energy. 7 These statements are independent of 
a power-law approximation. 

7 It is implicit in this paper that tt + and 7r~ responses are essentially the same. This 
was checked in a few cases, but was uninteresting. Since test beam work invariably 
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Fig. 7. Expected ir/p for one of the SDC reconfigurable absorber test-beam configu- 
rations [28] in which h/e « 0.67. to, Eq w , and Eq p are from the fits shown in Fig. 11 
of Paper I. 

In the power-law context, the energy-independent ratio f®-/fp should be 
(E 07T - 1 Eqp) 1- " 1 . This means that m is the same for both pions and protons. The 
Paper I (Fig. 11) simulations were consistent with equality. The scale energy 
E 0w - was found to be about 1 GeV, with some change from material to mate- 
rial. For protons the Monte Carlo simulations yielded E p o pa 2.6 GeV. These 
were consistent with the expectation that Eq was the approximate multiple- 
pion threshold[2]. Thus /°-//° « (1.0/2.Q) 1 "" 1 = 0.83 for m = 0.815 and 
fo/fo = 0.87 for m = 0.85. 

If h/e 7^ 1, a calorimeter should give a different response for charged pions 
than for protons. In the usual case, where h/e < 1, pions give the larger 
response. The effect is illustrated in Fig. 7, where as an example we use h/e = 
0.67, obtained from CALOR simulations (Paper I, Table 1) for the "iron" 
configuration of the SDC test-beam calorimeter [28]. It is regrettable that there 
was not time to measure the effect there. 

Equation (6) may be rewritten for the pion and proton cases: 




7r/e=l-(l-/i/e)/°_ 

P/e = 1 - (1 - h/e)f p (9) 

Rearrangement gives us the energy-independent ratio of the energy-dependent 
mean hadronic fractions: 8 



uses tt~ beams because they are free of p contamination, charged pions are simply 
labeled ir~ . 

8 This ratio was not calculated in Refs. [13] or [14]. 
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Table 2 

Calculation of the hadronic fraction ratio f®_(E)/fp(E) as a function of incident 
hadron energy (last column) using Eq. (11) and data from Table 2 of Akchurin 
et al.[14]. (/°-//°) = 0.859 ± 0.004, with X 2 = 13.2. Most of x 2 is contributed by 
the first and third points. 



Energy 


Response [14] 


f 7j- / fp 


(GeV) 


p/e 


ir/e 




200 


0.562 ±0.013 


0.647 ± 0.001 


0.806 ± 0.024 


250 


0.580 ±0.010 


0.648 ± 0.001 


0.838 ± 0.020 


300 


0.590 ± 0.006 


0.658 ±0.001 


0.834 ±0.012 


325 


0.592 ± 0.006 


0.652 ±0.001 


0.853 ±0.013 


350 
375 


fi n 7 +0.00i 
U - DU ' -0.004 

n R1 i +0.001 
U.OlX_ 00 3 


0.659 ±0.001 
0.664 ± 0.001 


0.868 ± 0.007 
0.864 ±0.005 


1 - 
1 - 


n/e 
p/e 







#-/■/? = T-^7T ( 10 ) 
^{E^-/E Qp f- m (11) 

The factor (1 — h/e) cancels. The constant ratios of energy-dependent quan- 
tities given by Eq. 11 does not depend upon a power law or any other model 
for the hadronic fractions, although the statistical sensitivity is maximal for 
small h/e. 

Specialization to the power-law case results in Eq. 11. Since only a = {l/h/e)El~ m 
can be found from ir/e (or p/e) measurements, the scale energy cannot be 
found without assumptions about h/e. It is therefore interesting that the ra- 
tio of the scale energies given by Eq. 11 can be very well known. 



Since Paper I was published, the CMS forward calorimeter group at CERN 
has measured the n/p ratio using a calorimeter consisting of quartz fibers 
embedded in a copper matrix (QFCAL) [13,14]. In this calorimeter, only light 
was detected, most of it coming from em showers, so that h/e was small, and 
the ir-p response difference was maximal, ir/e and p/e as a function of energy 
are reported in Table 2 of Ref. [14] and copied to our Table 2. 

This ratio is calculated in the right column of Table 2 for the Akchurin et al. 
data. These data together with estimates of the Paper I and a least squares 
fit to a constant /°-//~ are shown in Fig. 8. The fit yields f®-/f2 ~ 0.86. 

I am unable to find flaws in the several arguments leading to the conclusion 
that f®-(E)/ fp(E) is independent of energy, yet Fig. 8 shows evidence for 
energy dependence. A straight line with nonzero slope would certainly fit 
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Fig. 8. The mean hadronic fraction ratio f®-/f2 as calculated from the cop- 
per/quartz-fiber calorimeter data of Ref. [14]. The gray band is the range expected 
from Paper I; the crosshatched band is a constant value fitted to the data. The 
PDG scale factor is 1.6. 

the data better than a constant value. Akchurin et al. [14] needed to make a 
careful but difficult subtraction of n + contamination in their positive beam. 
The contamination was minimal at the highest energy, as is reflected in the 
uncertainties. If the pion contamination correction were overdone, one would 
obtain the observed low values at the lower energies. 

In the case of incident kaons, the leading hadron is probably a strange meson, 
but sometimes a pion. It is unlikely to be a proton or neutron. The response 
difference between incident pions and kaons should thus be small. 



5 mips 



As indicated above and in Fig. 1, e and h are the efficiencies with which 
electromagnetic and hadronic energy are converted into a visible signal. It is 
conventional to scale signal sizes, in ADC counts, to the mean response for 
minimum ionizing particles (mips), thus giving them something of an absolute 
meaning. In practice the "average" signal from penetrating muons, corrected 
for radiative losses, is assumed to be described by the Bethe-Bloch equation 
including the density effect. 9 It is then scaled to the value at minimum ion- 
ization, presumably defining the mip. 

But the mip is commonly used incorrectly. 

9 In everyday detectors, energy loss by escaping <5-rays or gain from entering 5- 
rays is small (2% level), so that "energy loss" and "energy deposit" can be used 
somewhat interchangeably. 
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In a simplification of the normal derivation of the Bethe-Bloch formula, 10 
ionization and excitation energy losses are calculated separately for (soft) dis- 
tant collisions (low energy transfer per interaction) and (hard) near collisions 
(high energy transfer). The regions are distinguished by the approximations 
appropriate to each [29-31]. One hopes for an energy at which they meet; this 
can sometimes be a problem in high-Z materials. Each contributes a factor 
In /?7 to the behavior at high energies: 

(1) As the particle becomes more relativistic, its electric field flattens and 
becomes more extended. This extension is limited by polarization of the 
material. This "density effect" asymptotically removes the ln/?7 factor 
contributed by the distant-collision region. The "relativistic rise" is still 
there, but with half the slope [32]. 11 

(2) The kinematic maximum energy T max which can be transferred in one 
collision sets the upper limit for hard energy transfer. Its rise with energy 
is responsible for the other In j3^f factor. As the particle energy increases 
there is more <5-ray production and the "Landau tail" grows and extends. 
The most probable energy loss, in a region well below minimum ionization 
that is dominated by many soft collisions, shows little or no relativistic 
rise and approaches a "Fermi plateau." This is more easily understood 
for the related restricted mean energy loss discussed below. 

The 6 rays are part of energy loss as described by the Bethe-Bloch equa- 
tion, and must not be confused with the density-effect correction 6(0^) 
or the muon radiative processes discussed below. 

For detectors of moderate thickness x (e.g., the scintillator tiles or LAr cells 
used in calorimeters), 12 the energy-loss probability distribution f (A; /3'j , x) 
is adequately described by the Landau (or Landau- Vavilov-Bichsel) disitribu- 
tion[33-35]. The most probable energy loss is 



ln^^ + lnf +J -^-^ 7 ) 



'121 



where £ = 0.153537 (Z/A) (x//3 2 ) MeV for a detector with a thickness x in 
g cm -2 , and j = 0.200. 13 While dE/dx is independent of thickness, A p /x 
scales as alnx + b. The density correction <5(/?7) was not included in Landau's 
or Vavilov's work, but it was later included by Bichsel[34]. It must be present 
for the reasons discussed in item (1) above. The high-energy behavior of S(P^y) 



10 Fano [29] introduces an intermediate energy transfer region. 

11 Review of Particle Physics 2006, hereafter RPP06. 

12 G < 0.05-0.1, where G is given by Rossi [[30], Eq. 2.7.0]. It is Vavilov's k[33]. 

3 Rossi [30], Talman[36], and others give somewhat different values for j. The most 
probable loss is not sensitive to its value. 
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Fig. 9. Bethe-Bloch dE/dx, two examples of restricted energy-loss rate, and the 
Landau most probable energy deposit per unit thickness in polystyrene scintillator, 
in which dE / dx\ m - m = 1.936 MeV g" 1 cm 2 . The incident particles are muons. 

is such that 



where hu p is the plasma energy in the material, 21.8 eV in the case of poly- 
styrene [RPP06, Tab. 27.1]. Thus the Landau most probable energy loss, like 
the restricted energy loss, reaches a Fermi plateau. The Bethe-Bloch dE/dx 14 " 
and Landau- Vavilov-Bichsel A p / x in polystyrene (scintillator) are shown as a 
function of muon energy in Fig. 9. It is interesting that the asymptote is 
nearly reached at 10 GeV for muons, and that it is not much higher than the 
minimum at just under 1 GeV. 

In the case of restricted energy loss [RPP06 Eq. (27.2)] the maximum kinetic 
energy transfer in a single collision is limited to some T cut < T max . One may 
find the energy- weighted integral of the <5-ray spectrum (d 2 N/dxdT; RPP06 
Eq. (27.5) ) between T cut and T max to find that it is equal to the difference 
between the restricted and Bethe-Bloch energy-loss rates. Similarly, one can 
integrate the 5-ray distribution over energy to find the number of 6 rays pro- 
duced in a tile — to find that in most CclSGS X dN/dx < 1. As the incident 
particle energy increases, the tail of the Landau distribution contains increas- 
ingly energetic but improbable energy transfers. Examples are shown in Fig. 9. 

In summary: The mean of the energy loss given by the Bethe-Bloch equation 
is ill-defined experimentally and is not useful for describing energy loss by 
single particles. 15 (It probably finds its best application in dosimetry, where 
only bulk deposit is of relevance.) It rises with energy because T max increases. 

14 1 follow convention and ignore the fact that dE/dx is actually negative. 

15 "The expression dE/dx should be abandoned; it is never relevant to the signals 



A 




2mc 2 ( 



(13) 



■p 



18 



1.0 




Electronic energy loss A (MeV) 

Fig. 10. Bichsel's calculations of the electronic energy deposit distribution for a 10 
GeV muon traversing a 1.7 mm silicon detector (which has roughly the stopping 
power of a 3-mm thick scintillator) [34,37-39]. The Landau- Vavilov function (dot- 
dashed) uses a Rutherford cross section without atomic binding corrections but with 
a maximum kinetic energy transfer limit T max . The solid curve was calculated using 
Bethe-Fano theory. Mo (A) and Mi (A) are the cumulative 0th and 1st moments of 
/(A), respectively. The fwhm of the Landau- Vavilov function is about 4£ for detec- 
tors of moderate thickness. A p is the most probable energy loss, and (A) divided 
by the thickness is the Bethe-Bloch mean, dE/dx. 

The large single-collision energy transfers that increasingly extend the long 
tail are rare. 

For a (3 fa 1 particle, for example, on average only one collision with T > 
1 keV will occur along a path length of 90 cm of Ar gas [3 9]. The energy- 
loss distribution for a 10 GeV muon traversing a 1.7 mm silicon detector, 
shown in Fig. 10, further illustrates the point. Here about 90% of the area 
(M (A)/Mo(oo)) but only ~ 60% of the energy deposition (Mi (A)/Mi(oo)) 
falls below the Bethe-Bloch (A), and at this energy A p fa 0.6 (A). The long 
tail of /(A) extends to T max = 4800 MeV. 

The mean of an experimental sample consisting of a few hundred events is sub- 
ject to these large fluctuations, sensitive to cuts, and sensitive to background. 
The mean (A) exp deduced from the data will almost certainly underestimate 
the true mean; in general A p < (A) cxp < (A). On the other hand, fits to 
the region around the stable peak in the pulse-height distribution provide a 
robust determination of the most probable energy loss. The peak is somewhat 
increased from the Landau A p by experimental resolution function. 

The ionization/excitation energy fraction sampled by the active region is con- 



in a particle-by-particle analysis." [39] 
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ventionally given by 

fraction sampled = t~j^ — • (14) 



S dE/dx 


scint 


SdE/dx 


scint + A dEj dx 


abs 



Boundary-crossing 5 rays are (accurately) assumed not to be important. Here 
S is the active region thickness, dE/dx | SC i nt is the energy loss rate in the active 
region (scintillator or other), A is the absorber thickness, and dE/dx]^ is the 
energy-loss rate in the absorber. It is "fairly sampled" if this is the case. It is 
not the case for electron and muon radiative losses. 

Relativistic muons also lose energy radiatively by direct pair production, 
bremsstrahlung, and photonuclear interactions. In iron at 1 TeV, these loss 
rates are in the ratios 0.58:0.39:0.03. The ratios are fairly insensitive to energy 
at energies where the radiation contribution is important. The pair:brems- 
strahlung ratio is about the same from material to material, while the pho- 
tonuclear fraction grows with atomic number. These contributions to dE/dx 
rise almost linearly with energy, becoming as important as ionization losses at 
some "muon critical energy" E^: 1183 GeV in plastic scintillator, 347 GeV in 
iron and 141 GeV in lead. 16 

The tables of Lohmann et al.[40] are commonly used. More extensive tables 
with a somewhat improved treatment of radiative losses are given by Groom 
et al.[31], and an extension to nearly 300 materials is available on the Particle 
Data Group web pages[41]. 1 

Muons in an absorber are accompanied by an entourage of photons and elec- 
trons (cascade products from direct pair production and bremsstrahlung) 
characteristic of radiative losses in the higher- Z absorber. If the calibration 
muon beam has been momentum-selected and then travels through air or 
vacuum to the calorimeter, it enters "naked," without its entourage of pairs 
and bremsstrahlung photons and, until this builds up over several radiation 
lengths, the signal distribution does not include the full radiation contribu- 
tion. If the equilibrium contribution is desired, absorber should be placed in 
the beam. In principle these radiative products should be detected with the 



Other charged particles experience radiative losses as well, but there is no easy 
mass scaling for the radiative loss rate. In a calorimeter incident high-energy pions 
lose energy by both radiation and ionization until they interact, but the higher loss 
rate is of little consequence and in any case the radiated energy is absorbed. 
17 For the PDG tables, an improvement to the pair-production cross section was 
made which slightly changes the muon dE/dx at high energies in high-Z materials. 
In both cases the ionization losses include a correction for muon bremsstrahlung on 
atomic electrons, so that at the highest energies the table entries slightly exceed the 
Bethc-Block values. 
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same efficiency as electrons, although a few high-energy pairs go through the 
active layers. There is little radiative loss in the low-Z active layers themselves. 



Bremsstrahlung is sufficiently continuous as to not introduce significant radia- 
tion fluctuations [42], but there are large fluctuations in the pair production en- 
ergy loss. Monte Carlo calculations by Striganov and collaborators [43] indicate 
that, while a high-energy shoulder appears on the energy-loss distribution, the 
most probable energy loss increases only slightly in "thin" absorbers," e. g. for 
1000 GeV muons incident on 100 g cm -2 of iron. They regard radiative effects 
as "important" when the most probable height of the normalized energy-loss 
distributions are lowered by <^ 10% when radiative effects are included. This 
is the case for the total signal from real calorimeters (more than 1000 g cm -2 ) 
at the highest muon calibration energies normally used. Although the most 
probable energy loss is still the best calibration metric, it does rise somewhat 
with beam energy because of the radiative effects. 



The calibration of the HELIOS modules (uranium/scintillator sandwiches) is 
particularly well described [44]. The common normalization of individual lay- 
ers was done via radioactive decay in the uranium plates, so the distributions 
shown are for entire modules. In correcting for the radiative losses they as- 
sumed fair sampling by the scintillator. Their correction of energy deposit for 
radiative effects is straightforward, but the robustness of the "average" energy 
deposit is unclear. 



Muon detection in SPACAL (scintillator fibers in a lead matrix) is carefully 
described by Acosta et al. [45]. The energy dependence of the most probable 
values is shown. The distributions clearly show both the radiative broadening 
and the increase of the most probable values due to radiative effects. 



In both of these cases, the calorimeters as a whole were calibrated. Since these 
are many radiation lengths thick, the lower average deposit at the beginning 
should not sensibly affect the result. 



All of this assumes muons of known energy. "Out of channel" muons, which 
have gone through or around the test-beam optics, are certainly not depend- 
able calibration particles, but are sometimes used [28]. Cosmic ray muons have 
a characteristic energy of about 3 GeV, but the flux falls off as about cos 2 9, 
where 9 is the zenith angle [RPP06, Sec. 24]. They can provide a useful if 
imprecise calibration in some situations. 
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6 e, h, and e/h 



Given a credible muon calibration, the quantity e/mip can be measured in 
an electron beam. In a sampling calorimeter, cascade electrons are predomi- 
nately produced and absorbed in the inactive higher- Z material, so the signal 
is significantly smaller than might be expected from the active layer's share 
of dE/dx (e/mip rs 0.6 to 0.7[46]), but with uncertainty associated with most 
probable energy deposit vs average energy deposit. It can be "tuned" by chang- 
ing the absorber /detector ratio, perhaps to achieve compensation (h/e = 1). 
Other things being equal, the detection efficiency is smaller if the absorber 
has a higher Z. The critical energy is lower, so characteristic shower electrons 
are more likely to deposit their energy before leaving the absorber. As a corol- 
lary, e/mip = 1 for a nonsegmented calorimeter (e.g., an inorganic crystal), 
and since there is always missing hadronic energy such a calorimeter is always 
noncompensating. 

The hadron efficiency h is more problematical; one finds a = (1 — h/e)E} ) ~ m or 
an equivalent by fitting the energy dependence of n/e, and assumptions must 
always be made about the constant multiplier to find h/e and hence h/mip. 
The multiplier Eq -171 is close to 1 for incident pions, but it is about 20% 
higher for protons, h is considerably more difficult to model, but in general it 
is smaller than e because of the wide variety of ways hadronic energy becomes 
invisible, e.g., through nuclear binding energy losses and "late" energy depo- 
sition (outside the electronics window). 18 It increases somewhat with Z, and 
can be enhanced by neutron production in uranium. It can also be "tuned" 
by the choice of material and sampling fraction [3]. 

Can we measure h/mip directly? Only by observing hadronic cascades in a 
calorimeter insensitive to 7r°-produced em cascades, or by observing cascades 
produced by hadrons below the 7r° threshold. In Paper I we speculated about 
building a calorimeter sensitive only to hadrons (a neutron detector) or to 
the em sector (a calorimeter sensitive only to Cherenkov radiation), but the 
context of the discussion was verification of the power-law approximation for 
and determination of the power m. 

In the spirit of only- (mostly-) em sensitivity in the copper/quartz- fiber CMS 
test calorimeter, Demianov et al. [47] made preliminary neutron measurements 
using Bonner spheres[48] adjacent to the copper/quartz- fiber calorimeter [13]. 
The longitudinal and transverse distributions were measured. The results were 
in fair to good agreement with MARS96[49] calculations, but not sufficiently 
detailed to obtain h/mip (or n/mip). Preliminary proposals [50] (in connection 
with International Linear Collider (ILC) detector R&D) are being made to 



A particularly nice discussion is given by Ferrari and Sala[17]. 
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measure the neutron flux by a variety of methods; future test-beam results 
will be of great interest. The problem will be discussed at more length in 
Sec. 9. 

One might use hadrons with energies below the 7r° threshold. ZEUS collabo- 
rators made measurements with low-energy protons and charged pions with 
a compensated U/scintillator calorimeter [51]. Interestingly, as the kinetic en- 
ergy of the beam was decreased from about 5 GeV to about 0.4 GeV, e/h 
decreased from its high-energy value (one) to the e/mip measured for elec- 
trons. The lower-energy particles tended to lose much or all of their energy by 
ionization, so they became indistinguishable from electrons at sufficiently low 
energies. The resolution also decreased from its hadronic value, approaching 
the em resolution until at the lowest energies noise became dominant. 

A more desirable (or complimentary) approach might be to use an incident 
beam of low-energy neutrons. Since Eq p ~ 2.6 GeV, one might expect the 
7T° threshold to be about T ~ 1.6 GeV. As the energy is scanned downward, 
a pure hadronic signal should emerge. The response would not be quite the 
hadronic signal observed from a higher-energy cascade, but this difference can 
probably be understood. At very least, measurements in a low-energy neutron 
beam would be interesting. The real problem is making the test beam. 



7 Resolution 

The arrows between boxes in Fig. 1 actually indicate the various p.d.f.'s de- 
scribing fluctuations in each of the steps. A more complete version, Fig. 11, 
defines these distributions, which are described in more detail in Table A.l. 
In the simple model considered in this paper, five p.d.f.'s appear. 

The potentially detectable energy deposit, or visible energy, is labeled "vis." 
It usually means ionization in a sensing medium, such as scintillator or liquid 
argon. In the rare cases where Cherenkov light is to be sampled, it means 
the Cherenkov radiation produced. It contains the variations due to energy 
deposit, not detection. The variance associated with the visible energy distri- 
bution at fixed E n o, dominated by fluctuations in the total kinetic energy of 
neutrons, is the intrinsic variance. 

This ionization is then sampled directly but more often via scintillators, where 
the scintillation light is usually detected by photomultipliers. The label "samp" 
refers to the additional fluctuations introduced in this process. 

The stochastic processes are defined as follows: 
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Fig. 11. Energy flow in a calorimeter, with the statistical distributions contributing 
to the experimental resolution indicated. 

(1) In the cascade initiated by the incident hadron, energy E w o is transferred 
to the 7T° sector via tt° production and decay. Because of its different 
energy dependence, em energy deposited by nuclear gamma rays is not 
included in my definition of E n o. The p.d.f. U(E n o) was introduced in 
Sec. 2 to describe its distribution. 

(2) The 7T° energy E v o is detectable with some average efficiency e via the 
ionization produced mostly by low-energy electrons and positrons. The 
p.d.f. of the em visible energy deposit at fixed E n o is labeled g^o(E^o\E 1T o). 

(3) Quite independently, the hadronic component produces a visible signal 
via energy loss by charged secondaries. This signal, the result of a variety 
of mechanisms, is produced with an overall average efficiency h. The vis- 
ible contribution by nuclear gamma rays is included here. Again, most of 
the ionization is by low-energy particles. The p.d.f. gh{E^ ls \E n o) describes 
the distribution at fixed E n o. 

(4) Only the total deposit E V1S = E^o+El 13 can be detected. Convolution over 
the intermediate variable E™ yields the p.d.f. -F vis (£' vls |£' 7r o). The prod- 
uct of this distribution with U(E w o), integrated over E w o, is the visible 
energy distribution which can be sampled. The width of this distribution 
is identified with the "intrinsic resolution." 

(5) Finally, the visible energy is sampled by measuring the ionization, ei- 
ther directly or by observing scintillation light with photomultipliers or 
photodiodes. 19 This step is to an extent under the control of the ex- 
perimenter, since it depends on scintillation efficiency, light collection 
efficiency, and other design details. For fixed -E V1S one measures a sig- 
nal E samp , chosen from a distribution F slimp (E slimp \E vls ) , which is then 
summed over the intermediate -E V1S to obtain the final distribution of 
the signal, E samp . Even this step is not a simple convolution, since the 
variance of F samp (S samp |£' vls ) is proportional to E V1S , not E. 

The intrinsic and sampling distributions were separated in a classic experi- 



The observation of Cherenkov light is another option. 



24 



ment by Drews et al.[52], who used compensating sandwich calorimeters with 
scintillator readout with either lead or uranium plates. Alternate sets of scin- 
tillators were read out separately. Sampling variations in the two sets were 
independent, while intrinsic fluctuations were correlated. These were recov- 
ered by adding and subtracting variances. 

Calculation of the combined distribution is tedious and not entirely obvious; 
the details are relegated to Appendix A. The result (repeating Eq. (A. 21)) is 

+ {l-h/efal {E). (15) 



Here cXg amp0 , (J 2 , and o\ Q scale the variances contributed by the sampling, 7r° 
energy deposit, and hadronic energy deposit, respectively. They have the units 
of energy. 

The first term is the familiar sampling contribution, except that it is multiplied 
by 7r/e. This is to be expected and required, since this contribution to the 
variance is proportional to the sampled visible signal, with mean (ir/e)E, 
rather than to the incident energy E. 

The two terms in the square brackets are the two pieces of the intrinsic vari- 
ance. Even if h/e = 1, the intrinsic variance has some energy dependence, 
since f n o increases with energy and fh decreases with energy. 

Together, the sampling term and the two intrinsic terms in the square brackets 
are usually represented as (C/VE) 2 , ignoring the energy dependence of each 
of the three terms. 

Wigmans[3] has noted that a- mtI /\^E for the simulated lead/LAr calorime- 
ter described in his Table 3.4 decreases with energy, reflecting the gradually 
increasing transfer of energy to the ir° sector. His calculated results for six 
energies, given in his Table 4.3, are plotted in Fig. 12. My curve was obtained 
by adjusting the intrinsic variance scales a 2 and cr^ . The best- fit parameters 
are a e0 = 5.1% and <Jho = 13.7%. The fit is remarkably good, and, as expected, 
aho is considerably larger than <7 e0 . 

The last term is the expected "constant term." Its mild energy dependence 
is discussed in Sec. 2. It must approach zero at high energies (1/E — > 0), as 
n(/jro) "crowds against" the f n o = 1 limit. That it can be represented as a 
constant is an artifact of the limited energy range of test-beam measurements. 
This point is discussed in more detail below. 

The 7T° part of the intrinsic variance fraction increases in importance as E 
increases, as does the sampling variance. Both curve downward if plotted vs 
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Fig. 12. Wigman's (simulated) data[3] for the intrinsic resolution of a Pb/LAr 
calorimeter scaled by l/y/~E, fitted with the term in the square brackets of Eq. (15). 
The fit parameter are given in the text. 
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Fig. 13. Shapes of contributions to sampling and intrinsic variance. As l/E — > 0, the 
slopes of the sampling and tt° contributions approach finite constants, since 7r/e — > 1 
and f^o — > 1, while the slope of the intrinsic hadronic contribution approaches zero 
(fh - 0). 

l/E, as shown in Fig. 13. The hadronic intrinsic contribution curves upward, 
since it decreases faster than l/E. 

It is difficult to verify Eq. (15), even with robust experimental data. The 
expected resolution should be a linear combination of the three curves shown 
in Fig. 13 (plus a constant term), so any deviation of the variance from the 
traditional Cf/E will show up as a slight curvature. Moreover, the sampling 
and 7T° contributions have such similar energy dependence that a simultaneous 
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Fig. 14. "Conventional" and Eq. (15) fits to the copper /quartz- fiber test module 
energy response data given in Table 3 of Ref. [13]. The fit parameters are given in 
Table 3. The fit shown by the dotted curve is for constant a w o, while the solid curve 
fit is for a more physically reasonable weak negative power law ^-dependence. The 
high-energy "turndown" is more clearly shown in (b), where the fits and data are 
shown relative to the usual straight-line fit. 

fit to cr samp0 and a e0 can be indeterminate. 

The square of the fractional energy resolution in the copper/quartz-fiber test 
calorimeter for incident pions (Akchurin et al.[13], Table 3) is plotted as a func- 
tion of 1/E in Fig. 14. The data curve downward relative to the "conventional" 
linear fit, Cf/E + Cf , shown by the dashed line. For this calorimeter intrinsic 
fluctuations were more important than sampling fluctuations except at the 
lowest energies, although sampling fluctuations were not negligible. Because 
of the nearly-degenerate shapes of the sampling and intrinsic ir° fluctuation 
curves, I set a samp0 = in making a fit, which is shown by the solid curve 
in the Figure. It describes the data well, and the physics responsible for the 
curvature is understood. 



Parameters for both cases are shown in Table 3. The fitted values for a e o and 
<Jho are very much larger than for the example discussed above and shown in 
Fig. 12; this follows from the excellent resolution of Wigman's model calorime- 
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Table 3 

Parameters for the three fits shown in Fig. 14. The power-law parameters m = 0.833, 
a = 0.753 and h/e = 0.247 (for Eq = 1-0) from Table 1 were used in the reduction. 
The experimental data are from the 3rd (a Tms /E) column of Tbl. 3 in Ref. [13]. 



a/E = d/y/E © C 2 Eq. (15) with cr samp o = Eq. (15) with cr sam po = 



Parameter 


Value 


Parameter 


Value 


Parameter 


Value 








377% 




372% 


Ci 


270% 




216% 




214% 






s±* 


14.2% 


Sl* 


15.7% 






S2* 


(fixed) 


s 2 * 


-0.058 




13.6% 


(1 - h/e) Sl 


10.7% 






X 2 /dof 


58.3/14 


X 2 /dof 


18.6/13 


X 2 /dof 


18.0/12 



* a n o = si (E/W0 GeV) S2 



ter and the (by design) poor resolution of the copper/quartz-fiber calorimeter. 

Other examples testing the 1/E dependence are hard to find. Many test-beam 
results are at low energy, many have large errors, and many of the earlier 
results are presented as functions of C\/\fE + C 2 rather than C\j VE © C 2 . It 
will be of interest to test Eq. (15) against further experimental results. 



8 Jets 

How is calorimeter's response to a jet different than its response to a single 
pion? There are three situations to consider: 

(1) An incident pion. The primary collision usually occurs about an inter- 
action length into the calorimeter. There is a minimum of backscatter 
("albedo"). The fragmentation process is dependent on energy and the 
nuclear environment. 

(2) A "test-beam jet," in which trigger counters ensure that the primary 
interaction occurs in a thin absorber in front of the calorimeter. This is 
exactly the same as the incident pion case, except for increased albedo 
because of the high probability that some of the first-collision debris 
interacts near the front of the calorimeter. 

(3) A primary fragmentation jet. The only evident differences from the above 
cases are the (much) higher energies and a simpler environment; except 
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in heavy-ion collisions just two particles interact. This section concerns 
whether the mix and distribution of photons, pions, and other particles 
results in calorimeter response different than the response to a pion or 
"test-beam jet." 

As elsewhere in this paper, the situation is highly idealized: The homogeneous 
or fine-sampling calorimeter is large enough to contain the entire cascade and 
the structure is uniform throughout. The realities of jet-finding and isolation 
algorithms, albedo, the effects of the magnetic field, passive material in front 
of the calorimeter, etc., are all ignored. 

The power law approximation for fh developed in Paper I will be used through- 
out this section. 

A jet with energy Ej consists of photons, mostly from 7r° decay, and "stable" 
hadrons. (Energy which might be carried away by leptons is ignored.) Since 
most of the incident "stable" hadron flux consists of charged pions, E « 
1 GeV and (1 - h/e) » a. 

One needs only to sum the calorimeter response to all of these particles to 
obtain the response to a jet. If R w 0j is the response to the jth 7r° (with energy 
E w 0j) in the jet and R hk the response to the kth. stable hadron (with energy 
E hk ), then the response to a jet is given by 

JV had 

X • X /''/,/, • (16) 

j=l k=l 



Using Eqns. (4)-(7) to evaluate Rhk and R n oj, this reduces to 



ET = eEj 



Vhad 

l-aEJ-^iEhu/Ej) 

k=l 



(17) 



Alternatively, the spectrum of stable hadrons can be described by a fragmen- 
tation function D(z), where z is the hadron's momentum parallel to the jet 
direction, scaled by the jet's momentum. In the present study z is treated as 
the fractional energy, i.e. z ~ E had /Ej. When the arguments leading to Eq. 17 
are repeated, one obtains 



eE, 



aE 



m—l 



l D{z)dz 



where D(z) describes the spectrum of all hadrons except for the 7r°'s. 
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The sum (Eq. 17) or integral (Eq. 18) thus appears as a correction factor 
to the normal hadronic response of a calorimeter. If it is unity, the response 
to a jet is the same as the response to a pion. In any case it is multiplied 
by a, which is usually < 0.3. The distinction between a single pion and a jet 
vanishes as the calorimeter becomes more compensated — except, of course, for 
the albedo, magnetic field, passive material in front of the calorimeter, and 
cone-cut effects mentioned above. 

If the sum or integral is evaluated for m = 0, the mean stable hadron multi- 
plicity (iVhad) is obtained. If m — 1, the result is the mean nonelectromagnetic 
fraction of the jet's energy (i^had)- The desired summation or integral, with 
m ~ 0.82-0.86, is in some sense an interpolation between the two. 

In using either experimental or Monte Carlo distributions to evaluate the 
sum or integral, special treatment of the very low- z region is necessary, as is 
normalization to an appropriate (-Fhad)- 

The integral in Eq. 18 is evaluated for four representative cases: 

Two experimental results, both with jet energies at or near Mz/2. Since the 
measurements are for charged hadrons, the distributions must be renormalized 
to include the contributions of such particles as A's and K^s. 

(1) Jets from Z decay, as measured by the DELPHI collaboration at LEP[53]. 
The published fragmentation function is for the entire event, so the func- 
tion has been normalized downward by a factor of two to describe the 
individual jets. Data were read from their Fig. 3(b) and extrapolated to 
z = 0. 

(2) CDF charged fragmentation function at y/s = 1800 GeV[54]. zdN ch /dz 
was extrapolated to z = to force (F,±) = 0.65, their reported value. 
(Since some of the energy is carried by neutrals, this value is probably 
too high for consistency with isospin conservation.) 

Two samples of TWOJET ISAJET[55] events at ^ = 40 TeV. 20 In both 
cases, all hadrons other than 7r°'s are used: 

(3) 3226 events with p t (hard scatter) > 40 GeV/c, and (100 GeV < Mjj < 
200 GeV). The mean jet momentum is 73 GeV/c, and the mean non-7r° 
hadronic multiplicity is 26. 

(4) 3042 events with p t (hard scatter) > 400 GeV/c, and 1000 GeV < Mjj < 
2000 GeV. The mean jet momentum is 677 GeV/c, and the mean non-7r° 
hadronic multiplicity is 70. The z distribution for these events is shown 
in Fig. 15. 



1 am indebted to my SDC collaborator E. M. Wang for running these simulations. 
This work was jointly reported in Refs. [6] and [7]. 
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Fig. 15. Distribution in z for ISAJET TWOJET events at yfs = 40 TeV, for 
1000 GeV < Mjj < 2000 GeV. 

The results are summarized in Table 4. There is ambiguity because of uncer- 
tainty in (.Fhad) m the simulations and (-F c h) i n the experimental results. If pion 
production dominates, one might expect (Fhad) ~ 2/3 from isospin consider- 
ations. (In Paper I, we reported fractions closer to 3/4.) Some of the bias can 
probably be removed by normalizing / zD(z)dz to 2/3, as indicated by the ta- 
ble entries in parentheses. As can be seen, the integral is slightly less than unity 
for the similar low-energy LEP and Tevatron fragmentation functions, and it 
is slightly greater than unity for simulated 40 TeV jets. Values lie between 0.84 
and 1.15 before normalization, and 0.92 to 1.06 after normalization — probably 
well within the uncertainty of the fragmentation functions in either the exper- 
imental or Monte Carlo cases. The integrals also change by about 0.05 if m 
is changed by 0.01, introducing an additional uncertainty which could be as 
great as 20%. Given the various uncertainties, I conclude that the correction 
factor for fragmentation jets at the highest-energy colliders should be between 
0.85 and 1.15. 

The compensation factor a ~ (1 — h/e) appearing in Eq. (17) and 18 serves to 
further reduce the effect of the correction factor in producing a jet/ir differ- 
ence. The percentage errors for the two limiting cases 0.85 and 1.15 are plotted 
in Fig. 16 for calorimeters with h/e = 0.70 and h/e = 0.85, values which might 
occur for a Pb/LAr or badly designed metal/scintillator calorimeter. The un- 
certainty in the exponent m could introduce an error of about 3% for jets 
below 100 GeV in a poor calorimeter. 

In the context of a power law approximation to the hadronic fraction for an 
incident pion, the response for an incident jet thus differs from the response to 
an incident pion by a simple correction factor, an integral over the fragmenta- 
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Table 4 

Integrals over representative fragmentation functions. Numbers in parentheses are 
calculated for the nonelectromagnetic energy fraction normalized to 0.67. In the 
case of the DELPHI and CDF results, the unrenormalized energy fraction is for 
charged hadrons only. 



Source 



Process 



DELPHI Z -> jet jet 11.0(12.1) 0.84(0.92) 0.61(0.67) 

CDF v / i=1.8TeV 17.8*(19.9) 0.94(0.97) 0.65(0.67) 

ISAJET 40 TeV, (pj) = 73GeV/c 26.2(25.2) 1.04(1.00) 0.69(0.67) 
ISAJET 40 TeV, {pj) = 677GeV/c 69.8(64.7) 1.15(1.06) 0.72(0.67) 
* The extrapolated low-momentum part of the function contributes 10 to this total. 
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Fig. 16. Energy determination error as a function of jet energy for representative 
values of h/e, for the two extreme case of the correction factor: 0.85 (top curves, 
for hard, low-multiplicity jets) and 1.15 (bottom curves, for soft, high-multiplicity 
jets). 

tion function. Given the uncertainties involved, no difference between jet and 
pion response can be found. 



9 Beating the devil 



An estimation of the ir° content of individual events would permit correction 
for intrinsic fluctuations (reduction of the "constant term"), along with its 
contribution to energy uncertainty. 

Several attempts have been made to use the radial and longitudinal detail to 
estimate, and correct for, the 7r°-induced cascades. During tests for the SDC 
construction, it was proposed that the 7r° contribution might come "early" 
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in the cascade, and could be estimated by excess energy deposit in the first 
layers. This turned out not to be true [56]. The ATLAS barrel calorimeter 
group adjusts downward the contribution of readout cells with large signals, 
since these tended to be from n cascades. In the test-beam runs they achieved 
slight improvements, e.g., from (46.9 ± 1.2)%/ y/E to (45.2 ± 2.2)%/ v / £[57]. 
Nearly a decade earlier, this algorithm was successful in correcting the WA 78 
data [20]. Ferrari and Sala have simulated such corrections for the LAr TPC 
ICARUS detector, where very detailed 3D imaging is possible, and conclude 
that the 7r° content can be fairly well determined from the direct observation 
of the em cascades [17]. The success of these corrections depends on the detail 
available, and any gains are usually marginal. 

Mockett[58] suggested long ago that information from a dual-readout calorime- 
ter with different h/e's in the two channels could be used to estimate the elec- 
tromagnetic fraction f n o for each event. Winn [59] has proposed using "orange" 
scintillator, observing the ionization contribution through an orange filter and 
observing the Cherenkov contribution through a blue filter. This has not yet 
been implemented, and looks problematical. 

The idea of using a quartz fiber/scintillator fiber dual readout calorimeter 
to extract an estimate of E T o for each event was discussed by Wigmans in 
1997[60]. Since then, the DREAM collaboration(Akchurin et al.[61]) has ele- 
gantly demonstrated the efficacy of the dual-readout technique, using a cop- 
per/optical fiber test-beam calorimeter. It consists of copper tubes, each con- 
taining three plastic scintillator fibers and four undoped fibers which produce 
only Cherenkov light. These are read out separately for each event. 

The principle is illustrated in Fig. 17. Akchurin et al.'s notation is used: S 
for the scintillator signal and Q for the Cherenkov signal, with both energy 
scales calibrated with electrons. For this example their values h/e\s = 0.7 
and h/e\Q = 0.2 are used for the scintillator and nominal Cherenkov readouts 
respectively. If there were only resolution contributions from the E n o distri- 
bution, events with different f^o should lie along the solid line drawn from 
f h = 1 to f h = (Uo = to Uo = 1): 

S = E(fo + fh{h/e\ s )) , and Q = E(f n o + f h (h/e\ Q )) . (19) 

The effects of finite resolution are illustrated using simulations which give 
f n o for 100 GeV negative pions axially incident on a very large lead cylin- 
der. (Results at 30 GeV from the same study are shown in Fig. 3(a).) For 
this cartoon example I arbitrarily introduced a Gaussian scatter in both E n o 
Ko/VlOO GeV = 1.5%) and E h (a h /VW0 GeV = 3.0%). The "events" are 
shown by the small +'s in the figure. The solid histogram at the top shows 
the marginal distribution in S. The mean is 84.7 GeV, the fractional standard 
deviation is 5.3%, and there is the usual skewness toward high energies. 
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Fig. 17. A toy model showing energy correction for 100 GeV pions in an idealized 
DREAM calorimeter, where Q is the response in the quartz-fiber readout and S is 
the response in the scintillator- fiber readout. The observed "events" are shown by 
the +'s, and the corrected events by the o's. Rotating to a frame in which the Q' 
axis is parallel to the event locus provides an equivalent reduction. 

Energy correction is straightforward. With the definition 

R= 1 -^^, (20) 
l-h/e\ s K } 



Eq. (19) can be solved to obtain 
RS-Q 

J-'COTT 7^ '-. • \ ) 

rt—l 



The circles in Fig. 17 show the same events after reduction via Eq. (21), and 
the dashed histogram shows the marginal distribution. The mean is 100.1 GeV, 
the fractional standard deviation is 3.4%, and there is no evident skewness. 
Complete compensation has been achieved by using the simultaneous readouts. 

Alternatively, a coordinate rotation to axes (Q 1 , S') can be made, such that 
the new Q' axis is parallel to the event locus. The projection of the event 
distribution onto the new 5" axis is of minimal width. Scaled upward by the 



31 



geometrical factor, it becomes the corrected distribution given by Eq. (21). 



The importance of a large "compensation asymmetry" is evident. If the stan- 
dard deviation in Q is oq and the standard deviation in S is as (both in GeV), 
then the variance in E corr is approximately 



R 



2 



Since i? can be well determined either from test-beam measurements of ir/e 
as a function of energy or from fits to the slope in a plot of Q vs S at one 
energy, the error in R has been neglected in writing Eq. (22). In the present 
example R = (1 - 0.2) (1 - 0.7) = 2.66, so o% = 3.20 erf + 0.36 a^. Given the 
Cherenkov readout, <Tq is likely to be much larger than erf. The price of the 
correction is an increased error on each event, but it is clear from Fig. 17 that 
there is compensating improvement. 

Alternate schemes simpler than DREAM would be desirable. Winn's scheme [59], 
taking advantage of the different colors of Cherenkov and (red) scintillation 
light, uses common detectors but still needs the doubled number of photomul- 
tipliers. Clean separation of the two signals would likely be difficult. LSND[62] 
used a weak scintillator and distinguished between directional Cherenkov light 
and isotropic scintillation light. But this was a very different kind of detector, 
a homogeneous low- Z detector used in a search for rare signals involving sin- 
gle, low-energy electrons. The electrons at the end of a high-energy shower do 
not remember the original direction. 

Descendants of DREAM are being studied under the rubric of the "Fourth 
Concept Detector" for the International Linear Collider [63]. One starting point 
might be the dual-readout quartz/scintillator DREAM concept. It can imme- 
diately be improved to get rid of signal correlation between adjacent fibers 
and simplified, e.g., by alternating scintillator and Cherenkov fibers in groves 
in sheets of absorber (possibly tungsten) which can then be combined as a 
sandwich. 

An interesting departure from the fiber calorimeter idea was presented by Zhao 
at a March, 2006 ILC workshop [64]. A "conventional" iron/scintillator-plate 
sandwich calorimeter is constructed in which lead glass tiles are substituted for 
all or part of the iron plates. The Cherenkov light is detected via waveshifter 
fibers in grooves in the tiles. Since lead glass (heavy flint glass) with specific 
gravities up to about 5.7 are available, the calorimeter thickness might not be 
that excessive. However, one expects that transverse structure observation is 
essential, and there remains the problem of neutron detection. 

Perhaps a better approach is to observe the fast, blue, directional, polarized 
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Cherenkov light from an inorganic scintillator. This is the object of present test 
beam work, using the scintillator PbWC>4[27]. (The slow component has only 
a 50 ns decay time, A max = 560 nm (yellow), but the scintillation efficiency is 
only 0.1% that of Nal.) 

But if neutron detection can be added to neo-DREAM, then the in-principle 
ultimate 1% hadron calorimeter resolution [60] might be approached. Ways to 
do this are under active investigation by the Forth Concept collaborators and 
others, who hope to exploit one or more distinguishing features of the neutron 
signal: 

(1) Neutrons distribute further from the core of a cascade than do other 
components. 

(2) Gamma rays from nuclear de-excitation following thermal neutron cap- 
ture are slow, in the several hundred ns range [3]. Given fast gating re- 
quirements, their signal is probably not useful. 

(3) In a hydrogenous scintillator, ionization from the proton recoils in n-p 
elastic scattering can be observed. 

(4) Most neutron detectors in nuclear physics take advantage of the large 
cross sections for the 10 B(n, a), 6 Li(ri,a), and 3 B.e(n,p) reactions, with 
the boron reaction being the most popular [65]. The gas-filled detectors in 
common use are impractical for a calorimetric applications, but boron- 
and lithium-loaded scintillators exist and are being further developed. 
The inorganic scintillator Lil(Eu) is an obvious candidate, but its crys- 
talline structure and 300 ns decay time both present problems. Organic 
borate additives in conventional plastic scintillators might have promise. 
There are common high-boron glasses and there are glass scintillators; 
whether a boron glass can be made to scintillate remains to be seen. But 
in any case, the time scale involved here is probably much too long. 

(5) Neutron interaction products are slow protons and fission fragments, so 
the nonlinear light output in scintillators (Birks' law) offers another, if 
unlikely, avenue. 

One might imagine planes of PbW04 functioning as a dual readout calorime- 
ter, with interleaved organic scintillator sheets. The PbW0 4 produces scintil- 
lation light from several ionization processes; the hydrogenous scintillator does 
the same (weighted a little differently), but it also detects the ionization from 
the n-p elastically scattered protons. For each event the PbW04 scintillation 
signal, the PbW04 Cherenkov signal, and the organic scintillator signal might 
be represented as a point in a data cube analogous to the two-dimensional 
"data square" shown in Fig. 17. It will be interesting to see if a correction 
formula as simple as Eq. 21 can be found. 

With sufficient work and a little luck, it seems likely that the Holy Grail of 
ultimate resolution will be implemented in future calorimeters. 
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10 Discussion 



The conceptual picture of the physics of a hadronic cascade and the scaling 
law it implies have been rich in consequences for understanding the behavior 
of a hadron calorimeter. The response ratio 7r/e is particularly simple, and the 
pion-proton response difference, in retrospect so obvious, was an unexpected 
surprise. If the incident hadron is a jet rather than a pion, the response is still 
given by Eq. (5), except that fh is multiplied by an integral over the frag- 
mentation function which appears to be near unity. The recent results using 
dual readout are explored and possible extensions are discussed. Resolution is 
described by considering in detail the various stochastic processes involved in 
a hadronic cascade; they cannot be simply convoluted. 
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A Resolution 

In the cascade initiated by a hadron with energy E, energy is transferred to 
the em sector via 7r° production and decay. The energy deposit in the resulting 
em cascades produce ionization with some efficiency e. Most of the ionization 
is via energy loss by the abundant low-energy electrons. In the case of the CMS 
developmental Cu/quartz-fiber test calorimeter [13], Cherenkov light samples 
part of the electron path length. Quite independently, the hadronic compo- 
nent produces ionization through the many mechanisms involved in hadronic 
cascades, again mostly by ionization by low-energy particles, with overall effi- 
ciency h. Each goes its stochastic way independently of the other. One must 
calculate the distribution of the sum of the contributions to the ionization 
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with the constraint that f w o is fixed, then integrate over f n o. Finally, the dis- 
tribution is "sampled" by directly collection ions or detecting scintillation light 
(or Cherenkov light) via photomultipliers or photodiodes. The resulting p.d.f. 
has a variance somewhat different than the usual a/E = a 1 '/y/E © cr const . The 
differences, discussed in Sec. 7, are easily understood physically 

To expedite the calculations, it is useful to associate a characteristic function 
(c.f.) (e mx ) = 4>{u) with each p.d.f. f(x). It is essentially the Fourier transform 
of the p.d.f., and is discussed in the Probability section of RPP06[32] and many 
other places [66]. Among the properties I will use are: 

• Convolution of p.d.f. 's becomes multiplication of c.f.'s: 

f(x) = J fi(x)f 2 (x - y)dy cj>(u) = (A.l) 

• Let the conditional p.d.f. of f<2(x\z) be 4>2(u\z) and the p.d.f. of z be fi(z). 
Then 

<f>(u) = J fx(z)(f) 2 (u\z)dz . (A.2) 

• If 02 (above) is of the form (j>2(u\z) = A(u) exp(ig(u)z), then 

<f>(u) = A(u)Mg(u)) . (A.3) 

where <f>i(u) is the c.f. of fi(z). 

• The c.f. of a Gaussian p.d.f. with mean m and variance a 2 is 

(j){u) = exp (imu - a 2 u 2 /2) . (A.4) 

• Higher moments may be included by continuing the series: 

4>(u) = exp (imu - a 2 u 2 /2 - i/j 3 u 3 /3\) + ... (A. 5) 

Here /X3 is the third moment of the distribution about the mean. The di- 
mensionless "coefficient of skewness" 71 = fi^/cr 3 was introduced in Sec. 2 
and will be used here. 

The arrows between boxes in Fig. (1) actually indicate the various p.d.f. 's 
describing fluctuations in each of the steps. A more complete version, Fig. 11, 
defines these distributions, which, along with their c.f.'s, means, and variances 
are given in Table A.l. The notation is somewhat verbose in the interest of 
clarity. Throughout the calculations, the primary hadron energy E is implicit 
and constant. Whether energies (such as E n o) or energies scaled by the incident 
energy (such as f n o = E n o/E) are used as variables is arbitrary. I make the 
split choice of using f n o and fh, and energies elsewhere, to prevent even more 
complex notation. 
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Table A.l 

Probability distribution functions (p.d.f.'s) and characteristic functions (c.f.'s) used 
in the resolution discussion. The primary energy E is an implicit conditional vari- 
able, and the em energy E w o is the independent conditional variable used in the 
development. The p.d.f. of the final sampled energy is not used explicitly. 



Distribution 


p.d.f. 


c.f. 


Mean 


Variance 


Fractional energy of 7r°'s 


nc/» 


0n(it) 


/tt° 




Ionization in em showers 
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Total ionization, fixed E w o 


F vis (S vis |^o) 
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Sampled signal, fixed E V1S 
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The p.d.f. n(/ 7r o) is discussed in Sec. 2. For reasons discussed there, f n o is 
chosen as the independent variable rather than its hadronic counterpart fh = 
1 — fir°- Typical simulations are shown in Figs. 2 and 3. The mean of f n o was 
defined as /° , the fractional variance was found to be a 2 , and its coefficient 
of skewness 71 was found to be about 0.6. Its c.f. is thus 



0n(«) = exp (m/° - «V 2 /2 - ra 3 7i^o3! + ...). (A.6) 



The skewness is carried forward in the calculation. The other p.d.f.'s are as- 
sumed to be near-Gaussian, with c.f.'s of the form given in Eq. (A. 4). 

The conditional p.d.f. gw (E™\fn ) describes the visible signal produced by 
the deposit of E f n o in the em sector. "Visible" means energy deposit, usually 
ionization, which can be sampled by an appropriate transducer. Its mean value 
is eE n o, 21 and its c.f. is (j) gn o(u\E n o). The variance for an ensemble of events 
with the same E w o should be proportional to E v o. The c.f. may be written as 

(j) gn o(u\E n o) = exp(iueE w o — u 2 al eE n o/2) . (A. 7) 

where a 2 scales the variance. Since the variance has units of (energy) 2 and is 
proportional to the energy, cr 2 has the units of energy. 

Similarly, the distribution of ionizing hadronic energy at fixed E n o is given by 
g^E^lE^o). The mean is h(E — E n o). (Since Eh = E — E n o, it is sufficient to 
express the condition on Eh as a condition on E n o.) In analogy to Eq. (A. 7), 
I write the c.f. as 

<f) gh (u\E n o) = exp(iuh(E - E n o) - u 2 a 2 h0 h{E - E^)/2) . (A.8) 

21 The normalization of e and h is ignored because at the end only the ratio h/e 
appears. 
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The complicated hadronic response is dominated by a small number of colli- 
sions with large nuclear binding energy losses, so its distribution is wider than 
the em response[24,44,67]. It is thus expected that ct 2 ^ > ct 2 , but as shown 
in Sec. 7 it is hard to distinguish the ways the contributions of ct 2 and the 
sampling term modify the energy dependence of the resolution. 

I interpret g^o^E^flE^o) and gh{EJ^\E^o) as the ir° and hadronic contributions, 
respectively, to the intrinsic resolution. This point will be explored later. 

Only the total ionization (or Cherenkov light) E V1S = E™ + E™ can be sam- 
pled. Let the conditional p.d.f. of E vls be F^E^lE^o): 

F vis (E vis \E n o) = J g n o(Ej\E n o) g h (E™ - E v j\E n o)dEj (A.9) 



This integral is a simple convolution, so by Eq. (A.l) 

(f) vis (u\E n o) = (f) gn o(u\E n0 )(f) gh (u\E n o) . (A. 10) 

The sum over E n o results in the distribution 

F vis (E vls ) = J n(/,o)F vis (r s |E ir o)^o . (A.ll) 

Via Eq. (A.2) the c.f. of F vis is 

0vis(«) = J n(/ w o ) vis (w \E w o)dE w o . (A. 12) 

The c.f. v i S (it | .E^o) can be calculated using Eqs. (A. 7) and (A. 8). For simplicity 
here and in the algebra leading to Eq. (A. 15), it is convenient to define Act 2 = 
CT 2 /i/e — CTgQ. Terms involving E^o are collected into the second exponential: 

vis (u|Ko) = e iuhE - u2 °™ hE/2 x e ^oMi-Ve)-» 2 A<x 2 /2) ( Al3 ) 
Written in this way, (j) vis (u\E n o) is of the form A(u) exp(ig(u)z), so by Eq. (A. 3), 

vis ( M ) = e iuhE-u*al hE/2 x n ( Me(1 _ h / e ) _ m 2 eA(T 2 /2) ? (A. 14) 

where g(u) is identified with e(u(l — h/e) — iu 2 Aa 2 /2). <pu_{u) is given by 
Eq. (A. 6), so it remains to substitute this function into Eq. (A. 14) and collect 
the terms multiplying powers of u. These terms can then be identified as 
the mean, variance, and skewness of F vis (E vls ). After considerable algebra, 
Eq. (A. 14) yields 
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u 



exp (me£(/° + fth/e) 
^eE[f n oa! + f h al Q h/e + (1 - h/efa^E 2 ] 
> 3 [ 7l ^ e 3 £ 3 (l - Ve) 3 + 3^ Aa 2 e 2 £? 2 (l - h/e) 



(A.15) 



The final step is to "sample" the ionization with whatever output transducer 
is being used. Although the experimenter has little control over the variance 
of F V i S (E vls ). 22 the design might be changed to improve light collection, for 
example, if the variance contribution due to photoelectron statistics were sig- 
nificant. Again a Gaussian distribution is assumed. The variance contribution 
from the sampling transducer is proportional to E V1S : 



W£ samp |£ vis ) 



9<7rrr 2 P vis 



exp 



(E 



samp 



E 



vis\2 



o 2 rpvis 
z ' u sampO^ 



'samp 



'samp 



U 



exp[zE™(u + f uV s 2 amp0 )] 
F^(E*\E 1t0 )<f> aam> (u\Er ia )dEr' 



(A.16) 

(A.17) 
(A.18) 



Since the variance is not a constant, a simple convolution is again insufficient. 



Following the recipe of Eq. (A. 3), g(u) 
in Eq. (A.15): 



u + ra 2 a 2 amp0 /2 is substituted for 



'samp 



[u + i«V 2 am] 



/2) 



(A.19) 



The mean pion response (the multiplier of iu in Eq. (A.15) is unaffected: 



V = eE(f° + f° h h/e) ; 
or vr/e = 1 - (1 - h/e) ft , (A.20) 

so that the usual form for n/e (Eq. (5) is recovered. 

However, eE (/° + f^h/ ^e)cx 2 amp0 is added to the variance of F wls (E wls ) (the mul- 
tiplier of —iu 2 /2 in Eq. (A.15)). The final fractional variance for the calorime- 
ter is 



a 1 



Cgo + A°Ve)a 2 amp0 
E 



E 



+ 



E 



h/efalo 



22 There are two caveats here: The effects of noncompensation can be minimized by 
the methods used by the DREAM collaboration[61], as discussed in Sec. 9 and (in 
principle so far) by measuring the neutron flux on an cvcnt-by-cvcnt basis[21] in 
order to reduce the intrinsic resolution contribution of gi % {E^ a \E 1 fi). 
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+ {l-h/e) 2 o$>{E), (A.21) 

where I have followed convention and scaled the energy to electron calibration: 
eE — > E. The energy dependence of a 2 is made explicit in the last line. This 
is almost the usual form for the resolution {a / E = C\/\[E © C2), but with 
some important differences: 

(1) The first term, the sampling contribution, is scaled by the em response. 
Since it is the ionization which is sampled, this contribution is perforce 
proportional to ir/e. 

(2) The terms in square brackets are the em and hadronic contributions to 
the intrinsic variance. The shape and interpretation of these terms is 
discussed in Sec. 7. 

(3) The analysis reproduces the familiar "constant term," with variance con- 
tribution explicitly proportional to (1 — h/e) 2 . Its important energy de- 
pendence is discussed in Sec. (2). 

The third moment about the mean (^3) of the sampled distribution is the 
coefficient of — iu 3 /3\ in (p sa , mp (u): 

fi 3 = 7 i^o£ 3 (l - h/ef + 3^ Aa 2 £ 2 (l - h/e) + 3£a 2 is a s 2 amp0 , (A.22) 



where the energy is again scaled to the electron calibration: eE — > E. Here 
Ecr^ is is the variance of F vis (E vis ), the coefficient of —u 2 /2 in Eq. (A. 15). 

The first term is to be expected in any noncompensating calorimeter; it is 
just the skewness of H(f n o) "playing through" to the end. As discussed in 
Sec. 2, the dimensionless coefficient of skewness, 71, is about 0.06 for the model 
discussed there (ir~ on Pb, using an old version of FLUKA), and a T o = 12.5% 
at 100 GeV with some mild energy dependence. 'jia^oE 3 is the actual third 
moment about the mean of H(f n o). 

It is interesting that the visible energy deposition and sampling terms also 
contribute to the skewness. In the first case, this is because the variance of 
the visible energy at fixed E n o is proportional to E n o, and so at large E n o a 
wider distribution is contributed to F v i S (E vls ) than for low E v o — even though 
for a given E n o the distribution is (taken to be) Gaussian. 

For the same reason, sampling also contributes to the skewness. The first two 
contributions both vanish if h/e = 1, but the third term does not. Even in the 
case of a compensating calorimeter, we should not expect an exactly Gaussian 
distribution. 
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